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1. Introduction 

Herein, we are interested in the homogenization of energy Uke quantities in elec- 
tromagnetism, and more particularly in Maxwell's equations, without boundary 
conditions. We use the notion of H-measures, introduced by Gerard and Tartar 
[5], [23]. We prove a general propagation Theorem for H-measures associated to 
Maxwell's system. This result, combined with the localisation property, is then 
used to obtain more precise results on the behaviour of H-measures associated to 
this system. 

As known, an H-measure is a (possibly matrix of) Radon measures on the product 
space fix 5 where ft C M" is an open domain and S is the unit sphere in K". 
In order to apply Fourier transform, functions defined on the whole of R" should be 
considered and this can be achieved by extending them by zero outside the domain. 
For this reason, we consider Maxwell's system in the full space R'^, which means 
without boundary conditions. Let us mention that similar works already exist, see 
in particuliar [2]. However, in [2], computations are far from being complete. 
If one is interested in coupling with boundary conditions, the usual pseudodifferen- 
tial calculus behind the notion of H-measures is not sufficient, and one should use 
much more technical tools. In the case of semi-classical measures, this is now well 
known, see for instance [6], [8], [12]. 

In the context of H-measures, and in particuliar without a typical scale, then one 
should use tools similar to those developped in recents works, see for instance [1]. 
However, the results presented here will be important for the full Maxwell problem, 
with suitable boundary conditions. 
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Let ft be an open set of M We will consider Maxwell's system in the material Q 
with electric permeability e, conductivity a and magnetic susceptibility fj given by 

i) dtD^{x, t) + J^{x, t) = iotH^{x, t) + F^, 

ii) dtB^{x,t) = -votE^{x,t)+G^{x,t), 
Hi) diYB^{x,t) = 0, 

divD^{x,t) = g^{x,t), 



(LI) 



IV } 



where X E il and t G (0,T), , , , and arc the electric, magnetic, 
induced electric, current density and induced magnetic fields, respectively. 
Morever, (the charge density), and are given, and we have the three 
constitutive relations 

f 1) D^{x,t) = e{x)E^{x,t), 



2) J^{x,t) = a{x)E^{x,t), 

3) B^{x,t) = rj{x)H^{x,t), 



(L2) 



where e, a and fj are 3x3 matrix valued functions and £ is a typical lenght going 
to 0. 

We shall consider this system in the full space Q = M^, without boundary conditions. 
Since we are not taking into account the initial data, we will also assume that the 
time variable t belongs to R. 

We shall use the notion of H-measure to compute for instance energy quantities in 
the following cases : 

i) Constant coefHcient case: here, we assume that the electric permittivity e, 
conductivity a and magnetic siisceptibility f/ are 3x3 identity matrices, i.e. 



e = a = 7i= (/d)3x3 • 



(L3) 



ii) Non constant coefRcient-scalar case: in this case, we consider that the 
matrix e, a, fj are scalar 3x3 matrix valued smooth functions, i.e. 

e = e(/(i)3x3, 0- = (T(/rf)3x3, fl = v{Id)3x3 (1-4) 

bounded from below. 



where e, a, r] are smooth functions, given in 
We will also assume that 



ft^^pt^Qty^ / -(0,0,0) in [L' 



p3\6 



X 



and 



{E^,H^y — in L^(IR x weakly. 



^) weakly (L5) 



(L6) 



After some prerequisites on H-measures, see [6] and [23] , presented in Section I, we 
use this notion in Section II to prove 
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Theorem 1.1. Constant coefficient case 

Assume (1.2), (1-5), (1-6) and (1-3). Then, up to a suitable extraction, the H- 
measure TTL = TTL(t,x, C); C = {Co,C'), C = {(1,(2,(3), associated to {u^), can he 
expressed as follows 

( (' ®('ait,x,() (' ^('c{t,x,()\ 

^= , , , , • (1-7) 

\C ®(d{t,x,() ( ®(h{t,x,()) 

Here a{t,x,() andb{t,x,() are positive measures, while c{t,x,() andd(t,x,() are 
complex measures such that c = d, all supported in [{Co = 0}U{C' = 0}]n{CiC2C3 = 
0}. 

They satisfy the transport system (propagation property) 

I C' P (-§ - 2a) -j2di'-{Tr{(' ® C'S' E)c] = 2Re Tr^^f,,, 
l'=i 

3 , „ 

J2 di'.[Tr{( ® ('d' E)a]- |C' l' ^ = 2i?e Tr^n/,,, 



(1.8) 



dt 

i'=i 

i C P (-^) + Y.^i'.\Tr{,( ® C'S' E)ci] = 2Re Tri^uf,,, 

v=\ 

- ^ d,,\Tr{( ® E)6]+ I C' P - 2d) = me TrM^/,,. 

Above, a derivative with an upper (resp. lower) index denotes a derivative wrt. 

variable (' (resp. x). E = E(C') is the constant 3x3 matrix whose action is given 
by E.a = (' A a, for all a G M!^. Finally, Huf is the 6x6 matrix correlating the 
sequences and f^, written with blocks of size 3x3. In particuliar, it is zero if 
at least f^ is strongly convergent to 0. 
Finally, one has the following constraint 

(idxia + (Idx^^a + (Id^^a = 2Re Tr^ugn 

and similarly for b, c and d, where ^i^g is the 6x6 matrix correlating the sequence 
with the sequence (P = (g^, 0, 0, 0, 0, 0)*. 

Theorem 1.2. Non constant coefficient- scalar case 

Assume (1-2), (1.5), (1.6) and (1.4)- Let the dispersion matrix (see formula (3.35) 
3 

below) he L{x,() = ^'^^ Aq^ {x)(jA-' , which has three eigenvalues, each with fixed 

multiplicity two, for (' ^ and given by 

Wo = ,uj+=v\( \ ,UJ- = -v\( |. 
Then the matrix P{x, () = Ao(CoId + L{x, ()) has also the following three eigenval- 



i^o = (o ,'^+ = Co +^'|C I ,1^- =(o-v\( 
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each with fixed multiplicity two, where v{x) = — is the propagation speed. 

Using the propagation basis and the eigenvector basis introduced in (3.51) and 
(3.52), it follows that the H-measure TTL associated to a suitable subsequence of 
can be expressed as: 



where TTLj^- are 3x3 matrix valued measures. Furthermore, one has 

TTLii = i[(C' «) C)ao + ^(^^ ® z'^)a+ + ^{z^ ® z^)b+ + ^{z^ ® z^)a^ + ^{z^ O z^)b-] 

TTlia = f [(^^ ® z^)a+ - {z^ ® z^)b+ - {z^ ® z^)a_ + (^^ ® z^)b_] 

TTL21 = f [(-2^ O z^)a+ - {z^ O z'^)b+ - o z^)a- + {z^ O z'^)h+] 

^ ^22 = ^[(C ® CO&o + \{z''® z^)a+ + ^1)6+ + ^ z^)a- + ^{z^ ® z^)b-] . 

using notations given by (3.51). Above ao, b^, a± and b± are all scalar positive 
measures supported in the set [{Co = 0} U {Co = \C \}\^ {C1C2C3 = 0}- Finally, 
one has the following propagation type system 

3 3 
-e{x)dt\[l^^ + Co ^ dve{x)d^'\[l^^ - 2c7TTln - a''E.5;,TTli2 = 2J?e/i„/^^ 
;'=i 

3 3 



-^(a;)Tni2 + Co 5^ a,.r?(ar)5''Tn,i2 + ^ ^'Ea^TTln = 2iieM„/,, 

i'=i 

3 3 
-£(x)at 11121 + Co ^ di,e{x)d''T\l^^ - 2aTn2i - ^ 5''E.5,.Tn22 = 2i?eM„/,, 



i'=l v=i 

3 3 



-r?(^)Tn22 + Co ^ 5i-r,(a;)9''m22 + ^'Ea^-THai = 2i?eM„/,, 
i'=i i'=i 

w/iere we are using same notations as in Theorem 1.1 for the right hand side. 



2. Some basic facts on H-measures 

In this Section, we recall some results from the H-measures theory, taking the 
presentation of Tartar [23]. However, this is also similar to the exposition of Gerard 
[5,7], relying upon Hormander [10], [11]. 

Definition 1. Let il bo an open set of M " and let bo a sequence of functions 
defined in M " with values in M^'. We assume that converges weakly to zero in 
(L^(M"))^'. Then after extracting a subsequence (still denoted by e), there exists 
a family of complex-valued Radon measures (TTL C))i<i.j<p on M" x S'^~^, 
such that for every functions 4'i,4'2 in Co(M"), the space of continuous functions 
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converging to zero at infinity, and for every function in (^(S* " the space of 
continuous functions on the unit sphere S"~^ in M", one has 



<TII,^-,0i</>2®V'>= / / <PihHC/\C\)~^ijix,OdxdC 

= hm / [Fid^^umwihu^mm/mc 

e > jRn 

Above, F denotes the Fourier transform operator defined in i^(R"), for an inte- 
grable function / as [F{f)]{() = / f{x) e-'^'^'^-^dx, while F is the inverse Fourier 

JR" 

transform defined as [F{f)]{x) = / f{() e^'^^'^''^ dC,. z denotes the complex conju- 

Jr" 

gate of the complex number z. 

The matrix valued measure TTL = (TTL,ij)i<i<p is called the H- measure associated 
with the extracted subsequence . 

Remark 1. H-measures are hermitian and non- negative matrices in the following 
sense 



TTUj = TTLjj and 

P _ (2.2) 

~^3't'i't>o > for all <j)i (j)2 (t>u e Co(M"), 

and it is clear that the H-measure for a strongly convergent sequence is zero. 

Although we consider the scalar case for all properties of H-measures, all the fol- 
lowing facts are easily extended to the vectorial case. 

Definition 2. Let a e (7(5 ""i), 6 G Co(M"). We associate with a the linear 
continuous operator A on L-^(M") defined by 

FiAum = a{C/\C\)F{um a.e. C € (2.3) 

and with b we associate the operator 

Bu{x) = b{x)u{x) a.e. a; e M" . (2.4) 

A continuous function P on R" x 5""^ with values in R is called an admissible 
symbol if it can be written as 

P{x, = Y1 ^»(^) ^ = J2 bnix)an{C) (2.5) 

n=l n=l 

where a„ arc continuous functions on 5""^^ and 6„ are continuous bounded func- 
tions converging to zero at infinity on R " with 

max|a„(C)| max |6„(a;)| < oo . (2.6) 
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An operator L with symbol P is defined by: 

1) L is linear continuous on i^(M") . 

2) P is an admissible symbol with a decomposition (2.5), satisfying (2.6). 

3) L can written as the following form 

n 

L = AnBn + compact operator 

n=l 

where An , Bn are the operators associated with a„ , 6„ as in (2.3) and (2.4). 
With notations as in (2.3) and (2.4), one can show that the operator C := AB — BA 
is a compact operator from L^(R") into itself. Denote by X"*(IR") the space of 
functions v with derivatives up to order ra belonging to the image by the Fourier 
transform of the space L^(R") i.e. (_F(L^(R"))), equipped with the norm 



Then, if A and B are operators with symbols a and b as in (2.3) and (2.4), satisfying 

one the following conditions 

1) a e C^{S "-!) and b G X\W) , 

2) a e Xl^{W \ {0}) and b G CqHM"), 

it follows that the operator C = AB — BA is a continuous operator from L^(IR ") 
into i?^(R") and extending a to be homogeneous of degree zero on R", then 

VC = (AB - BA) has the symbol 

OXi 

(VrV.,C^C.i:^^. (2.7) 

K=l 

The main results of H-measures theory are given by the next two results 

Theorem 2.1. Localisation property Let be a sequence converging weakly to 
zero in (_L^(R"))^' and let TTL be the H-measure associated to . Assume that one 
has the balance relation 



n 

Y^—{A^u^)^0 {Hr^l{n)r strongly, 
fe=i 

where A'' are continuous matrix valued functions on Q C M Then, on fl x S , 
one has 

n 

P{x,OT(l={J2CkA''{x))T(l = 0. (2.8) 
fe=i 

This result shows that the support of the //-measure TTL is contained in the (char- 
acteristic) set 

{{x, QenxS''-^ , det P{x, = 0}. 



I 



I 



1, 2008 21:3 WSPC/INSTRUCTION FILE H-measures 



Theorem 2.2. Propagation property for symmetric systems Let be given 
matrix valued functions A'^ in the class Cq (O) . Assume that the pair of sequences 
{u^, f^) satisfies the symmetric system 



fe=l 



(2.9) 



and that both sequences (u^), (/^) converge weakly to zero in L'^{Qy. Then the 
H-measure p, associated to the sequence {u^,f^) and given under the form 



Mil Ml2 

M21 A'22 



(2.10) 



satisfies the equation 

n 

< /xii, {P, ijj+i^Yl ^''^'^ - >=< '^^iTrpi2), ^ > (2.11) 



fe=i 



for all smooth functions tp{x, (). Here S := 1/2(B + B*) is the hermitian part of 
the matrix B and {P,tlj} is the Poisson bracket of P and tp, i.e. 



{PM = d^Pd,^ -d^^d.p . ±i§^ - 1^1^). 

^ dC,i dxi dQ dxi 



(2.12) 



3. Applications to Maxwell's system 

This section is devoted to the proofs of our main results stated in the Introduction. 



3.1. Proof of Theorem 1.1: Constant coefficient case 

This case corresponds to the assumption (1.3), that is all the matrices e ^i) ,and a 
are the identity matrix, i.e 



e = a = rj 



'1 0' 
1 
.00 1 



[Id) 



3x3 



In this case, system (1.1) can be rewritten as 



i) 



dt 



-{x,t) + E^{x,t) = TotH^{x,t)+F^{x,t), 



a) 

Hi) 
iv) 



{x,t) 



dt 

Hi) d\vH^{x,t) 
dxvE^{x,t) 



= -votE^{x,t)-\-G^{x,t), 
0, 

p^{x,t). 



(3.1) 



(3.2) 
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Recalling the notation of the Introduction, it follows that Maxwell's system (3.2) 
can be written as 



and 



Here 



and 



t^-'i^-f- (») 



■^°=(:"l = f"") (3.5) 



fy y V Id 



The constant antisymmetric matrices Qk , 1 < < 3 are given by 
/CO 

Qi = -1 I , Q2 = I I , Qa = I 1 I (3.7) 
\0 1 

the matrix C and by 

'a 0\ _ /Id 0\ fF^ 

oo) [ 

Matrices i?% i = 1,2,3, are given by 











where the 3x3 matrices /3* are such that 

/3ii = except for = 1 
Finally we have denoted = {g^, 0, 0, 0, 0, 0)*. 

Denote the H-measure corresponding to (a subsequence of) the sequence by 

(3.9) 

The measure TTL is a 2 x 2 block matrix measure, each block being of size 3x3. 
In the following, let xq = t, x = {xo,x), x = {xi,X2,X3). We let C, denote the dual 
variable to x, with C = (Co, C), C = (Ci, C2, Cs)- 

To state the localisation property (2.1), we need first to express the symbol of the 
differential operator appeating in (3.3), for which one has 



I 



I 
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and thus 



P{x, C) = CoA'ix) + CiA\x) + C2A\x) + C3A^{x) = 



-Id o\ / qn / Qi 

Idi+^MQi J+^nQ2 



, . , Ql\ fCold -E 

+43 



where 



Qa ; V E Cold 



E= Cs -Ci . (3.10) 




Clearly E is antisymmetric (i.e. E* = — E), so that P is a symmetric matrix. 
Using the localisation property, it follows 

'Cold -E 

PTTl= I 11 1=0 (3.11) 



E Cold 




me ' rri 



and thus 



f 1) CoId.//e+Et.i/„,e = 0, 

2) CoId.l^em+Et.Z/^ = 0, 

3) E.I^e + Cold.i^em = 0, 
. 4) E.Z/em + Cold.fm = 0. 



(3.12) 



First note that from (3.11), since (see also next subsection) P has Co) ± | C' I 
eigenvalues, that TTL is supported in {Co = 0} U {C' = 0}. Then, one has 

Lemma 1. The H-measure TTL can be written under the form 

( C' (E)C'a{t,x,0 C' (^C'c{t,x,0\ 

\ , , , , (3-13) 

\C (E>Cd{t,xX)C (E>Cb{t,xX)J 

where a, b are scalar positive measures, c and d are scalar complex measures such 
that c = d, all supported in {Co = 0} U {C' — 0}. 

Proof of Lemma 1 
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Multiplying (3.12-1) par Co, (3.12-3) par E* and substracting the results, one has 

(CoId + E2).i/e = 0. (3.14) 
We discuss the following distinct cases 

i) case Co = 0. Then note that C since C belongs to the unit sphere of M^. 
From (3.12), we have E.i/e = 0. 
Then, we use the following Lemma 

Lemma 2. 

// E.A = 0, then the matrix A has the form ^ = C ® for some vector a G M ^. 
Proof of Lemma 2 

We denote the columns of the matrix A by the vectors 

A = [ai di ai] . (3.15) 

But as E.A = 0, we get that 

[Eai Eai Ea{\ = (3.16) 

or 

Eai = ,i= 1,2,3. (3.17) 

For i = 1, and similarly for the other cases, ai = ( /'2 | G M^. Then from (3.10), 
(3.17), we get that 

-Cs C2 \ / Vl\ /-W2C3+f3C2 



Cs -Ci U2 = ^^iCs-^^aCi |=C«)ai=0 (3.18) 

-C2 Cl Q / \V3/ \-WlC2+W2Cl, 

which implies that Q / /a\ and thus a\ = c\C, , for some constant c\ G M. 
Thus all in all, and for i = 1,2, 3, all the columns of the matrix A are parallel to 
the vector C = (Ci) C2) Cs); so we can write that 3^ = CiC, , and by arranging these 
numbers Cj as components of the vector a, we get that 

A = a^C'. (3.19) 
End of the proof of Lemma 1 

Using Lemma 2, we can conclude that I'e = C ® C c-it^ 2;, C ) and thus the blocks of 
the matrix if-measure, which satisfy system (3.12), are such that 

t'e = C' O C'a{t,x,Q, 
t'm = C'OC'H*,a;,C), 

(3.20) 

fem = C' C'c{t,X,C), 
t'me = C' ® ('d{t, X, C) = i^em = C' ® C'c{t, X, C) 



I 



I 
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where a{t, x, () and b{t, x, () are real positive measures, while d{t, x, () and c{t, x, () 
are scalar complex measures, such that c = d. 

ii) case Co ^ but C ' = (Ci, C2, Cs) = 0. 

From (3.14) one has ((^Id.iy^ = (recall that ( =0 implies that E = 0), and in 
this case we have Cc^-Ve = 0. As (o ^ 0, wc get — 0. Repeating the same steps, 
for the other equations of the system (3.12), we get that f e = t'm = ^'em = t'me = 
and finally, one has 




(3.21) 



iii) case Co 7^ and C' = (Ci> C2) Cs) 7^ 0- Using (3.14), one has 

(Co^Id + E.E).!/e = <^ Q.v^ + E.(E.z/e) = 0. (3.22) 

But as (E.t'e) = C' <8i t'e, one has that E.(C' ® v^) = ® (C' ® and from (3.22), 
one gets 

Co^-Z^e + C ® (C ® J^e) = 0. (3.23) 

Now we shall show that fe = 0. 

By contradiction, we assume that ^ 0. Using (3.23), one has (1 — C ^)-z^e = 
-C <H) (C ® i^e), thus either C'^ = 1 or C'^ < 1 or C'^ > 1- 

If C ^ = 1, one has Co = 0, which is a contradiction with Co 7^ 0. If C ^ > 1, then one 
have Co^ = 1 — C ^ < which is not possible, because Co^ > 0, and if C ^ < 1 then 
we get again Cq^ = C ^ ~ 1 < 0) which is a contradiction. 

Thus all in all, we conclude that = ^ and similarly for v^n = J^em = t^me = 0, 
and finally we have also, in this case 




ending the proof of Lemma 3.11. 

Let us now turn to the localisation property associated with equation (3.4). In this 
case, it follows that TTL, satisfies 

3 

(X^B^O)Tn = o. 

3 

Setting S(C) = ^ B^Q, since det B{Q = C1C2C3) it follows that TTL is supported 
J = l 

in the set {C1C2C3 = 0}. 
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All in all, the scalar measures a, b, c and d are all (also) supported in the set 
{C1C2C3 = 0}. 

Now, wc wish to write down the propagation property, and for this purpose, we 
need to compute the Poisson bracket, associated with equations (3.3) and (3.4). 
Letting ip = 'tjj{x,Q he a,n arbitrary smooth function, recall first that the Poisson 
bracket is given by 



{P, V'} = X! - d^ijdiP 



_y^ dPd-ip dip dP 



(3.24) 



Recall that a derivative with an upper (resp. lower) index denotes a derivative wrt. 
variable ( (resp. x). 

In our case, starting with (3.3), we have 



( ^ ' \ 

i'=i 

3 

C^d''E.di'ip) Iddtip 

\ i'=i 



(3.25) 



Next, we compute the term -0 ^ diA'' - 2ipS, where S = 1/2(C + C *) = C. Note 

also that as A' , < Z < 3 , are constant matrices, one has ipdiA'' = 0. 
Adding all the terms, we obtain 



{P,V} + ^V'5;A'-2V'S: 



;=o 



Iddttp -(^a''E.9i'V) 

i'=i 

3 



i'=i 



fipld0\ 


V / 



and thus 



{P,i>} + J2^diA^-'^i^S: 



I ^ ' \ 

3 

(J2d^'^.di'tp) Iddttp 



(3.26) 



I 



I 
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Taking into account the form of TTL deduced from the locaHsation property (3.11), 
we get 



3 U ®Qa{t,xX) C OCc(t,a;,C) 

< TTL, {-P, ^} + ^ ^diA^ - 'i^S >=< 

1=0 \C' <E)C'd{t,x,C) C <^Cb{t,x,C) 



V 



> 



(3. 



{dti) - 2^)id -(^ ^A'^) ^ (^-^^^ 

i'=i 

3 

i'=i 

so that the propagation property reads as 

3 

< C' ® C'ait, X, C){dt^p - 2V')Id + C' «) C'c(i, x, C)(^ 5''e.9;/ V) >= ^Re^iuf,, , 

/'=i 

3 

< -C' O C'a(i,a;,C)(X] '^''e.S^'V) + C' O C'c(t, x, C))Id5tV >= 2ReHuf,^ , 
i'=i 

3 

< C' «) C'<^(t, a;, C)(5t^ - 2V')Id + C' C'K*, a;, C)(5I ^''e.S;'^) >= 2i?eM„/,, , 

;'=i 

3 

i'=i 

with the notations explained in the statement of Theorem 1.1. 

Writing these equations in V and then taking the trace of each equation, we get 

finally (1.8). 

Now, it remains to take into account the propagation property coming from (3.4). 

3 

Since in this case, the characteristic polynomial is given by B{() = ^ -B-' Cj, it 

J = l 

follows with a small computation that one has the following propagation property 

/ C' <8 C'a{t, X, c' ® C'c{t, x,0\ ( r(x, r(5;, c) \ 
< , > = '2.Re^lue 

VC ®Qd{t,x,0 C ®Qb{t,x,Q)j \T{x,QT{x,Ol 

again using the notations in the statement of Theorem 1.1. Here T{x,C,) is the 3x3 
matrix given by 

r(i,C)= diag{dxi'4',dx2'4',dx3tp). 

Writing each equation and taking the trace, this gives the last constraint men- 
tionned in the statement of Theorem 1.1. 
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3.2. Proof of Theorem 1.2: Non constant coefficient- scalar case 

Let us recall that we assume (1.4), that is e, i) and a are 3x3 scalar matrix valued 
functions given by 



(e{x) \ 



e = e(Id)3x3 = 



e(x) 
e{x) 



(3.29) 



and 



V = '7(Id)3x3 



(ri{x) \ 



V 



ri{x) 
r]{x) 



,'& = CT(Id)3x3 



((7{x) \ 

cr(.X') 

a{x) 



(3.30) 



where e, r] and a are smooth functions in C^(M^), bounded from below. 

The first equation of Maxwell's system can then again be written in the form of a 

symmetric system 



where 



AO 



x)^ + C{x)u^ = f^ 



e{x)ld 
T]{x)ld 



(3.31) 



(3.32) 



is a (2 X 2) block matrix with 3x3 blocks, the constant antisymmetric Q k ■ 1 < ^ < 3 
and the matrix A' , 1 < i < 3 being the same as in (3.6) and C now given by 



C{x) 



a{x)ld 




(3.33) 



is a (2 X 2) block matrix, each block being of size 3x3. 

As the sequence converges weakly to zero in L'^{Q.)^, again up to a subsequence, 

it defines an H- measure TTL, which is a 6 x 6 matrix valued measure. 

As in the preceding constant case, to express the localisation property linked with 

(3.31), we compute the associated symbol P(x, Q which is here given by 



Coe(x)Id -E 
E Co??(a;)Id^ 



(3.34) 



with E is still given by (3.10). The localisation property then states that PTTL = 0. 
Let us first show the following 



I 



I 
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Lemma 3. Assume that C' 7^ and let L{x, Q = ^^^0 ^i^)Cj'^'' be the dispersion 
matrix. Then L has three eigenvalues, each with constant multiplicity two, given by 

Wo = 0, uj+ = Co + v\C\, w_=Co-w|C|. 
The matrix P'{x,Q = (^old + L{x,(^') has also three eigenvalues, given by 

wo = Co = Co + v\C'\ = Co -■J^IC'l, 

1 



each with constant multiplicity two, where v{x) 
speed. 



\Mxjr]ix) 



is the propagation 



Proof of Lemma 3 

L{x,( ) can be rewritten as 



L(x,c') = EV(^)o^^ = - 

3=1 



or in block form 



Lix,C) 



/ 




Cs/m -C2//X 

-C3/M Ci/ii 
V C2/M -Ci/m 



-1/eE 




-Cs/e C2/e \ 

<i/e 





Cs/e 
-C2/e Ci/e 










/ 



(3.35) 



The action of the matrix E = E(C') is also given as 

E(C')[p]=C'Ap 



(3.36) 



(3.37) 



where E(C ) is given in (3.10), for all p G M.^. Letting iv be an eigenvalue of L{x, C ) 
corresponding to an eigenvector X, one has 



-1/eE 



or 







l//iE 



-1/eE 




(3.38) 



(3.39) 



Let us check that w = is an eigenvalue. In order to see this, we have to solve the 
algebric system 



-l/eE(C )v = 0, 

i/mE(c')u = o 



(3.40) 
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which is equivalent to 

f C A V = 0, 

(3.41) 

It follows that u and v arc colincar to C', and thus that (0, C) and ((' , 0) is a basis 
of the eigenspace corresponding to the eigenvalue of L. It shows also that is 
indeed an eigenvalue of multiplicity two. 

It remains to find the eigenvalues u ^ of L. For this purpose, set 

We can as well assume that <J G S"^, the unit sphere in R^. Then a; = =f1 are 

eigenvalues of the matrix D(( ). 

Indeed, we have to solve the system of equation 

— A V = LOU , 

(3.42) 

C' A u = wv . 

For a; = 1 to be an eigenvalue of the matrix D{( ), we have to solve the system 

-C' A V = u , 

(3.43) 

C' A u = V 

which admits a 2D space of solutions, and for a; = —1 to be an eigenvalue of the 
matrix D{( ), we have to solve 

-C A V = -u , 

(3.44) 

which again admits a 2D space of solutions. Thus, all in all, we have obtained that 

a;i=— l,Ci;3 = l 

are the eigenvalues of multiplicity m = 2 of the matrix D{( ). Recall that co 2 = 
is also an eigenvalue of the matrix £)(C ) with multiplicity m = 2. 
Next, note that if A is an eigenvalue of a matrix A then cA is an eigenvalue of the 
matrix cA, where c > is a constant . Thus we conclude that 

-|C'l,0,|C'| (3.45) 
are the eigenvalues of the matrix B{( ) = L{x, C ) with multiplicity m = 2, for all 

Set ^{x, C) = A(a;, C ) ~ Co and recall that if X{x, () is an eigenvalue of the matrix 

P corresponding to an eigenvector u = u{x,(^), then /i(.T, C) = A(x,C ) — ^0 is an 
eigenvalue of the matrix L[x,(^ ) corresponding to an eigenvector X = X(a;, C ). 




I 



I 
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If /x(a;, Q = X{x, C) — Co is an eigenvalue of the matrix L(x, ( ) corrsponding to an 
eigenvector X = X(x, (J ), then A(x, () = n{x, () + Co is an eigenvalue of the matrix 
P'{x,C) corresponding to an eigenvector X = X(x,C)- As the eigenvalues of the 
matrix L{x,C, ) are given in (3.45), we can conclude that 

Co-«IC'l ,Co ,Co + «|C'l (3.46) 

are the eigenvalues of the matrix P'{x,() with multiplicity m = 2. Now we shall 
show that the propagation speed v is given by 

y(^X^ ^ 

Indeed, from (3.39) and for w ^0, one has 

f-l/eE(C')v = '^u 
I l/lJ,E{C')u = u-v 

or 



(3.47) 




(3.48) 



and thus 



Thus 



— C' A u . 



C'A(^C' Au) = -ec.;u: 
C' A (C' A u) = -e^w^u : 
{(' A iC A u)}.u 
(t Au).(C Au) = -e/xw2||-||2 
IICAulp 

but C' is orthogonal on the vector u and thus 



e/iw^||u|P 



-e/xa;^||u|P 



(3.49) 



,2 iiciniui 



UJ = ± 



IICII 



± 



IC 



(3.50) 



This ends the proof of Lemma 3. 

For the eigenvectors of the matrix P'{x, (), one first chooses an orthonormal basis 
of R*^. This basis consists of the propagation triple in the direction of propagation 

C and of the two transverse unit vectors z^^''{('),z^^''{('). 

Let (C) -z^^^HO) -2^'^HC')) ^ be this basis of propagation. In polar coordinates 
they are given by, see for more details [15], [16] 

' sin 6 cos 4)\ I cos Q cos < 



c 




.(1) 



(C) 



COS f sm ( 
— sin< 



.(2) 



(C) 




(3.51) 



I 



I 
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Where Id = (C? + Cl + Cl)'/'- 

Then, one can show that the eigenvectors of the matrix P'{x, () are given by 

The eigenvectors and 6q represent the non-propagating longitudinal and the 
other eigenvectors correspond to transverse modes of propagation with respect the 
propagation speed v. 

Lemma 4. The H-measure TTL, has the form 

m, = 6o (g) boQo + 6o ^ blbo + b\(» b\_a+ + 6+ (g) 6+6+ + 6l_ 6La_ + 6?_ (g) 6?.6_(3.53) 

where oq ,6o two positives measures supported in the set {(^q = 0}, a+, 6+ are 
two positives measures supported in the set {C,q = —v \ C,' |}, and a_, 6_ are two 
positives measures supported in the set {Qo = +v \ C,' \} . Here b^, b^, b\, b\, b^_, 
and 6?_ are the eigenvectors of the matrix P' given by (3.51) and (3.52). 

Proof of Lemma (4) Since PTTL = by the localization property, it follows that 

(^o)"'m = 0. (3.54) 

Since P'{x,C) = {Aq)-^P{x,C,), thus (3.54) becomes 

P'(a;,C)TTL = 0. (3.55) 

It follows that the support of the H-measure TTL is included in the set 

?7 = {(a;, C) e R'' X 5^ , detP' = 0} . (3.56) 

Note that, for every t = (^o,x,Q fixed, the matrix P' is diagonalizable. In fact, we 
will discuss the following cases : 

a) Case of equal eigenvalues 

i) If w = w + , (resp.w o = w _), then from Lemma 3), one has C,' = 0, (resp.w o = 
uj _), and thus Co = ±1, and P' = (old, which is diagonal. 

ii) If Wo 7^w+, to ^ to but w + = w _. In this case, again from Lemma 3, one 
has C = 0, and thus Co = ±1 and P' = (old, which is diagonal. 

b) Case of distinct eigenvalue 

Ifwo ^ CO - , then again from Lemma (3), using the basis of eigenvectors 

corresponding to an the eigenvalues wq, a; +, a; _, given by (3.51) and (3.52), P' = 
Cold is diagonal. 



I 



I 
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Next, as the support of the H- measure TTL is contained in the set of points (3.56), 
and as the matrix P' is diagonahzable, recalhng that the determinant of the matrix 
P' is given by (eventually with powers) 

detP' wow+W- . 

It follows that the support of the H-measure TTL is included in the set 

U = {{x,C) eR'^ X ,uo = 0[Ju+ = OUw.} . 

Next, one has 

P'TTl,. =0 j = l,2, 6 



where TTLj is the j — th column vector of the matrix TTL. Using Lemma (3) , we get 
that, for Wo = 0, there exist two scalars aj , fij such that 

Repeating the samcs steps for the other eigenvalues ui ^ , - , and using the hermi- 
tian property of the H-mcasurc, this ends the proof of Lemma 4. 

In order to write down the propagation property for equation (3.31), we first need 
to compute the Poisson bracket, which in this case, is given by 



{P,i^} = 



i'=i 

3 , 3 , 



.(3.57) 



3 



It follows that, since dkA^ = 



k=0 



{P,i;} + i;J2^kA''-2i;S = 
fe=o 

3 ^ 3 , 

e{x)IddtiJ -CoJ2^ ^'^|J^l>e{x)Id - 2^o-Id ''e.^^ ' V) 

i'=i i'=i 

3 , 3 , 

(^ 9 E.a, > V) v{x)Id^t^|J - Co 5] 9 r?(ar)Id 



(3.58) 



i'=i 



Writing the H-measure TTL as 



/Ta.n TrLi2 



(3.59) 



I 



I 
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where Tl\^j are 3x3 matrix valued measures, it follows that one has 

3 3 

i'=i i'=i 
3 3 



''=1 3 (3.60) 



l'=i 

3 3 



-?7(a;)Tn22 + Co dMx)d''TJl22 + Yl ^''Ea^TTLai = 2i?eM„/,, 
i'=i i'=i 

which is the form given in Theorem 1.2. 

Using the eigenvector basis (3.51) and (3.52), with the decomposition given by 
Lemma 4, we note that the elements of the H-measure TTL defined in (3.59), can be 

expressed as 

TTlii = i[(C' C')ao + ^iz^^ z^)a+ + ^{z'^ O z^)b+ + ^iz^0 z^)a- 



mz^^z^)b-], 

TTl,i2 = f [(2^ «) z^)a+ - (^2 (g, z^)b+ - (zi (g) 22)a_ + (2^ O zi)6_] 
TTI21 = ^[(-z^ 2^)a+ - (zi (g) z^)b+ - (z2 (g zi)a_ + (z^ ® 
m,22 = ^[(C' ® CO&o + 1(^2 ® + i(zi zi)6. 



(3.61) 



^LV5 . gV-^ ^ 2 ^ / + 

+i(z2 «) z2)a_ + i(zi (g) z^)6_]. 

If one wants to find an equation for ao for instance, one can proceed as follows. 
Recalling the form of TTL;^^ just given above, we take the equation for it in (3.60) 
and apply it to the vector 

Finally, because the divergence constraint in Maxwell's system only involves scalar 
valued functions, similar statement as in the end of Theorem 1.1 holds true again. 
Acknowledgements: The author would like to thank Radjesvarane Alexandre 
for several discussions and suggestions during the preparation of this paper. 
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